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Аннотация 

Введение. Многокритериальную оптимизацию с учетом противоречащих друг другу критериев задействуют для 
улучшения эффективности производства, сокращения затрат, повышения качества продукции и экологической 
безопасности процессов. В литературе описано использование многокритериальной оптимизации для 
производственных целей, в том числе при выборе условий реакции и улучшении технологических процессов. В 
представленной работе объект исследования — это процесс гидрирования полициклических ароматических 
углеводородов (ПАУ) при получении высокоплотных топлив. Для определения оптимальных условий процесса 
решается задача многокритериальной оптимизации на основе кинетической модели. Критерии: максимизация выхода 
целевых нафтенов и конверсия исходного сырья. Цель работы — создание программы, реализующей алгоритм 
многокритериальной оптимизации М№ЗСА-П (англ. поп-дотлтае зогипе оепейс а]еогит П). Благодаря этому на 
основе кинетической модели можно рассчитать оптимальную температуру для процесса гидрирования ПАУ. 
Материалы и методы. Для решения многокритериальной задачи оптимизации применялся генетический 
алгоритм М5СА-П. Используется также измененный отбор родителей и выживания в рамках фронта Парето. При 
необходимости разделения фронта решения выбирались по манхэттенскому расстоянию между ними. Программа 
реализована на языке Руфоп. 

Результаты исследования. В системе обыкновенных нелинейных дифференциальных уравнений химической 
кинетики концентрацию обозначили у» условное время контакта реакционной смеси с катализатором — т. Систему 
решили для реакции гидрирования полициклических ароматических углеводородов. Расчеты показали, что при т = 0 
У(0) = 0,025; у2(0) = 0,9; уз (0) = 0,067; уэ(0) = 0,008; ук0) = 0, 1=3—5,7,8,10-20; О(0) = 1. В качестве управляемого 
параметра рассматривали температуру процесса по двум критериям оптимальности: максимизация выхода целевых 
нафтенов (Д) в конце реакции и максимизация конверсии исходного сырья (ф). Значения /\ были в границах 0,43-—0,79; 
конверсии — 0,01-0,03; температуры — 200-300 К. Рост температуры сопровождается увеличением выхода целевых 
нафтенов и снижением конверсии исходного сырья. Каждое полученное решение — неулучшаемое. При 
моделировании процесса гидрирования ПАУ запустили алгоритм с размером популяции — 100, количеством 
поколений — 100. Разработана программа, реализующая алгоритм М№ЗСА-П. Рассчитано оптимальное множество 


значений температуры реакции гидрирования ПАУ, позволяющее получить неулучшаемые значения критериев 
оптимальности — максимизации выхода целевых нафтенов и конверсии исходного сырья. 

Обсуждение и заключение. Алгоритм №СА-П эффективен для решения задачи недоминирования и вывода 
оптимального решения для всех критериев. Будущие исследования следует посвятить подбору оптимальных 
параметров алгоритма, позволяющих увеличить скорость решения. Опираясь на полученные теоретические 
оптимальные условия реакции гидрирования ПАУ, можно реализовать процесс в промышленности. 
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Введение. Оптимизация многостадийных реакций применяется в химической, нефтегазовой, пищевой и других 
отраслях промышленности. На практике задачи оптимизации — многокритериальные, причем критерии зачастую 
противоречивы и имеют оптимум в разных точках. Многокритериальная оптимизация сохраняет актуальность, так 
как позволяет учитывать несколько параметров и выбирать наилучшее решение из множества вариантов. 

В рамках данной работы объект исследования — каталитическая реакция гидрирования полициклических 
ароматических углеводородов (ПАУ). Они представляют собой класс органических соединений, молекулы 
которых содержат не менее двух бензольных колец [1]. ПАУ распространены в межзвездной среде, входят в 
состав тяжелых фракций нефти, образуются при лазерном облучении углеродных материалов. Изучение этих 
соединений интересно с точки зрения выявления зависимостей между их химическим строением и физико- 
химическими свойствами. Кроме того, данные, полученные в результате таких научных изысканий, можно 
задействовать для создания новых органических и гибридных соединений с прочным углеродным каркасом, 
которые применимы в наноархитектонике. 

Прикладная наука соотносит наличие ПАУ с целями производства. Например, желательно присутствие ПАУ 
в сырье, если оно используется для получения кокса с заданной структурой [2]. Однако при производстве топлива 
ПАУ может негативно сказаться на эксплуатационных характеристиках продукта, например на плотности [3]. 

К выработке высокоплотного реактивного топлива предъявляются крайне серьезные требования. При 
высокой плотности оно должно иметь температуру кипения не выше верхней границы температуры кипения 
керосиновой фракции. Другой обязательный критерий — низкое содержание ароматических углеводородов. 
Отметим также дороговизну известных технологий получения высокоплотных топлив. 

Учитывая сказанное, задачу необходимо решать по принципу доминирования Парето для определения 
множества неулучшаемых вариантов с применением генетического алгоритма недоминируемой сортировки [4]. 

Цель исследования — разработка программы, реализующей алгоритм многокритериальной оптимизации 
М$ЗОА-П (Моп-доттаеа зогпе вепейс а1еогИбт П) и позволяющей на основе кинетической модели рассчитать 
оптимальную температуру для процесса гидрирования ПАУ. 

Материалы и методы 

Математическая модель. Опишем изменения концентраций компонентов в зависимости от времени 
реакций. Для этого используются уравнения химической кинетики, которые представляют собой систему 


обыкновенных нелинейных дифференциальных уравнений (СОНДУ): 
/ 


Чу: р 
Вы ) Уи, =1,...[, (1) 
ат & 
= 
1 1 Е — 
(09 Ви и И 
м) -®] [со "| [со "К; = А; еф| -—— | Ку =А--ефр|--= |. (2) 
ВТ ВТ 
1 1=1 
Здесь начальные условия: ух0) = у? при т = 0; те [0, 1]; у; — концентрации реагентов, молярные доли; 
т — условное время контакта реакционной смеси с катализатором, кг‘мин/моль; / — количество стадий; 
[— количество веществ; у, — стехиометрическая матрица; у’; — скорость /-й стадии, 1/мин или моль/(кг-мин); 
К, К; — константы скоростей стадий (приведенные), 1/мин; а; — отрицательные элементы матрицы уу; 
В; — положительные элементы уз; А;, А; — предэкспоненциальные множители, 1/мин; ЕЁ’, Е; — энергии 


активации прямой и обратной реакций, ккал/моль; А — газовая постоянная, 2 кал/(моль’К); Т — температура, К; 
т’ — продолжительность реакции, кг-мин/моль. 

Модель каталитического гидрирования ПАУ учитывает динамику мольного состава и объема реакционной 
смеси. Следовательно, принимаются во внимание изменения концентрации компонентов в каждый момент 
времени [3]: 


49 хо 4» и 
ат я 90-0, 


1 
1 1 
м, ] =, -] [(25№. 
а 9 О 


Для описания нестационарной реакции, которая происходит с изменением объема реакционной смеси, 
необходимо решить систему нелинейных дифференциальных уравнений в каждый момент времени. Прямая 
кинетическая задача представляет собой решение СОНДУ (1)—(3). 

В процессе гидрирования полициклических ароматических углеводородов из исходных ароматических 
углеводородов получают нафтены, которые обладают более высокой плотностью и могут использоваться в 
качестве топлива для реактивных ракет. С этой целью применяются никелевые катализаторы, а управляющим 


(3) 
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или варьируемым параметром является температура процесса, которая должна быть в пределах 200-500 К. 
Критерии оптимальности — максимизация выхода целевых нафтенов в конце реакции и максимизация конверсии 
исходного сырья. 

Алгоритм исследования. Многокритериальная оптимизация — это выбор из множества альтернатив 
наилучшего решения с учетом нескольких критериев. Важность каждого из них определяется весом 
(приоритетом). 

Допустим, /(х) — это целевая функция, а ограничения, заданные в виде равенств й1(х)...й»(х) и неравенств 


8„-1(х)...8>(х), представлены вектором-столбцом компонент х = [х1,..., Хн|7 в и-мерном Евклидовом пространстве. 
Сформулируем задачу нелинейного программирования [5]. 
Оптимизировать 
У(х) > ежи, хе Е", (4) 
при т линейных или нелинейных ограничениях в виде равенств: 
й;(х) = 0, 1 =1,...,т, (5) 
и при (р — т) линейных или нелинейных ограничениях в виде неравенств: 
2:(х)> 0, /=т+1,...,р. (6) 


Декомпозиция (4)-—(6) представляет собой постановку и решение задачи линейного и квадратичного 
программирования. Каждая из них определяется видом уравнений (4)—(6). Так, в случае квадратичной функции (4) 
и линейных уравнений (5), (6) — это описанная ниже задача квадратичного программирования. 

Определим экстремум функции 


1(х) = ао +с'х+х' Ох —> ем, (7) 

с ограничениями: 
а'х>Ь,х>0. (8) 
В уравнениях (7), (8) О — неотрицательно определенная квадратичная симметричная матрица; 


а, Ь, с — матрицы коэффициентов. 
При постановке многокритериальной Паретовой задачи оптимизации (4)—(6) будет иметь вид: 


ехёЕ (х) = Е(х")= Е". (9) 
В уравнении (9) Ех =(л (х), 2 (Х),..., 1" (Х)) — вектор-функция критериев оптимальности Л и ХВ. 


Множество х* — искомое решение задачи в области параметров варьирования. Множество Е” представляет собой 
искомое решение задачи в области критериев оптимальности, неулучшаемое в смысле аппроксимации Парето. 
Тогда х" определяет множество Парето, Е” — фронт Парето. 

Для решения (9) применили априорные и апостериорные алгоритмы аппроксимации Парето. Один из 
них — метод идеальной точки, которая представляет собой лучшее решение по всем критериям [6]. Чтобы ее 
найти, нужно сначала определить минимальные и максимальные значения каждого критерия всех 
рассматриваемых решений. Затем для каждого критерия выбирается максимальное значение из всех 
минимальных и минимальное из всех максимальных. 

Однако у данного подхода есть недостатки. Во-первых, он может быть неэффективным, если идеальная точка 
находится вне области допустимых значений критериев. В таком случае используются другие методы решения 
задач многокритериальной оптимизации. Кроме того, метод идеальной точки не учитывает взаимосвязь между 
критериями и может привести к выбору компромиссного решения, которое не является оптимальным по всем 
критериям. Поэтому при использовании данного метода необходимо дополнительно анализировать и проверять 
оптимальность полученных решений [7]. 

В решении задач многокритериальной оптимизации используется также метод лексикографического 
упорядочивания. При этом критерии упорядочиваются по приоритету и рассматриваются последовательно. Если 
решения не могут быть отсортированы по первому критерию, то они сортируются по следующему критерию 
ит. д. [7]. Преимущества лексикографического упорядочивания — простота и прозрачность. С помощью этого 
метода можно получить единственное оптимальное решение, которое легко интерпретируется. Что касается 
недостатков, отметим, во-первых, невозможность учесть компромиссные решения, которые могут оказаться 
оптимальными по всем критериям. Во-вторых, есть риск выбора невыгодного решения, если первый критерий 
имеет большой вес, но не является самым важным для данной задачи [8]. 

В представленной работе приводится решение задачи многокритериальной оптимизации процесса 
гидрирования полициклических ароматических углеводородов с применением известного метода М$СА-П. Он 
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основан на генетическом алгоритме и использует несколько техник для решения проблемы недоминирования [9]. 
Ниже описаны основные шаги алгоритма [10]. 

1. Инициализация популяции. Начальная популяция создается случайным образом. 

2. Оценка популяции. Каждый элемент популяции оценивается по нескольким критериям. 

3. Сортировка популяции. Элементы популяции сортируются по уровню недоминирования. Доминирующие 
и не испытывающие доминирования (недоминируемые) элементы помещаются в первый уровень. Элементы, 
доминируемые только элементами первого уровня, помещаются во второй уровень, ит. д. 

4. Выбор родительских элементов. Для создания новой популяции выбираются родительские элементы из 
нескольких первых уровней. 

5. Кроссинговер и мутация. Родительские элементы проходят кроссинговер и мутацию, чтобы создать новые 
элементы популяции. 

6. Оценка новой популяции. Новые элементы оцениваются по критериям. 

7. Сортировка новой популяции. Новые элементы сортируются по уровню недоминирования. 

8. Выбор новой популяции. Из новой популяции выбираются элементы для следующего поколения. 

9. Повторение шагов 4—8 до достижения критерия остановки. 

М$СА-П позволяет точно и оперативно работать с задачами многокритериальной оптимизации. Он эффективно 
решает проблему недоминирования, что позволяет получать оптимальные решения для всех критериев. 

М5СА-П основан на генетическом алгоритме с отбором родителей и выживания. Особи выбираются по 
фронтам, при этом фронт разделяется, если не все особи могут выжить. Решения в разделенном фронте 
выбираются на основе расстояния между ними, которое является манхэттенским расстоянием в пространстве 
критериев [9]. Крайние точки сохраняются на каждом поколении и им присваивается условно бесконечное 
расстояние для использования в последующих итерациях [11] (рис. 1). 


Рис. 1. Визуализация фронта Парето и выбора решений на основе расстояний 


На рис. 1 — пример множества решений для задачи многокритериальной оптимизации с использованием 
критериев Л и №. Красные точки обозначают фронт Паретто. Показан расчет расстояния скопления для 
решения 7 — это средняя длина стороны кубоида, в котором находится решение 1 (отмечено синей рамкой). 

Для усиления воздействия на отбор родителей М5СА-П использует двоичный турнирный отбор [9]. Каждая 
особь сначала сравнивается по рангу, а затем по расстоянию между ними. 

Результаты исследования. Программа, реализующая алгоритм решения задачи многокритериальной 
оптимизации, написана на языке Ру®оп. 
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В ходе исследования решена система дифференциальных уравнений (1-3). Для реакции гидрирования 


полициклических ароматических углеводородов она имеет вид [12]: 
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Прит = 0; у1(0) = 0,025; у>(0) = 0,9; уб(0) = 0,067; уэ(0) = 0,008; ук0) = 0, [= 3—5,7,8,10-20; О(0) = 1. 
На рис. 2 приведено рассчитанное множество решений, удовлетворяющих ограничениям системы, и фронт 
Парето, где /! — это выход целевых нафтенов, а р — конверсия исходного сырья. 
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Рис. 2. Фронт Парето процесса гидрирования полициклических ароматических углеводородов 


Управляемый параметр — температура процесса [13]. Критерии оптимальности — максимизация выхода 
целевых нафтенов в конце реакции и максимизация конверсии исходного сырья [14]. Расчеты по предлагаемому 
алгоритму дали результаты соответствия критериев оптимальности и значений параметра температуры, которые 
обобщены в таблице 1. 

Таблица 1 
Соответствие критериев оптимальности и значения параметра температуры процесса гидрирования 
полициклических ароматических углеводородов 


Выход целевых нафтенов /1 Конверсия исходного сырья № Температура Т, К 
0,43 0,03 200,00 
0,59 0,03 221,20 
0,71 0,01 250,00 
0,76 0,01 271,10 
0,79 0,01 300,00 


Итак, при увеличении температуры выход целевых нафтенов увеличивается, а конверсия исходного сырья 
уменьшается. Каждое из этих решений является неулучшаемым, выбор конкретных значений зависит от лица, 
принимающего решение. 

Для успешной работы №МСА-П необходимо правильно выбирать параметры алгоритма. В частности, следует 
оптимизировать критерии размера популяции, количества поколений, вероятности кроссинговера и мутации. 
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При запуске алгоритма на модели процесса гидрирования полициклических ароматических углеводородов 
использовались следующие параметры: размер популяции — 100, количество поколений — 100. Недостаточно 
большая популяция может обусловить преждевременную сходимость алгоритма к локальному оптимуму. 
Слишком большая способна замедлить процесс оптимизации [15]. Чрезмерное количество поколений 
потенциально ведет к переобучению алгоритма, тогда как при недостаточном количестве поколений может не 
хватить времени для достижения оптимальных решений. 

Обсуждение и заключение. Создана программа, реализующая алгоритм многокритериальной оптимизации 
М5ОА-П. Работа с соответствующей задачей в рамках этого метода включает решение системы 
дифференциальных уравнений, визуализацию множества решений, удовлетворяющих ограничениям системы, и 
построение фронта Парето. Кроме того, найдены значения варьируемых параметров для достижения целей 
оптимизации. Для процесса гидрирования ПАУ на основе кинетической модели рассчитано множество значений 
температуры, оптимальных для получения неулучшаемых значений двух критериев оптимальности: выход 
нафтенов и конверсия исходного сырья. С повышением температуры растет скорость реакции и выход нафтенов. 
Однако уменьшается конверсия сырья. К тому же слишком высокие температуры могут стать причиной 
побочных реакций и разложения продуктов. 

Данные, полученные в рамках представленной работы, могут быть полезны для оптимизации процесса гидрирования 
ПАУ в промышленных условиях. Важно учитывать влияние температуры на выход нафтенов и конверсию сырья при 
разработке стратегии производства. Кроме того, следует принимать во внимание и другие параметры, от которых зависит 
кинетика реакции. Это, например, давление, скорость потока реагентов и роль катализаторов. 

Таким образом, разработанная программа и предложенный алгоритм позволяют проводить одновременный 
анализ нескольких критериев оптимальности процесса на основе кинетической модели и формировать множество 
неулучшаемых значений варьируемых параметров. 
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